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Abstract 

If produced at the LHC, long-lived charged particles (LLCPs) would leave tracks in the muon 
detector. Time-of-Flight based methods for detecting these particles become less efficient if the 
LLCPs are fast, which would typically be the case if they are produced in the decays of some 
mother particle which is either heavy or very boosted. Thus for example, in supersymmetric 
models with long-lived sleptons, the long-lived sleptons produced in neutralino decays are often 
fast, with P > 0.95 even at a 7 TeV LHC. We propose to use the (mis-measured) invariant mass 
distribution of "muon"-lepton pairs, where the "muon" could be a slepton LLCP, to detect it. 
This distribution peaks somewhat below the neutralino mass. The peak can be further enhanced 
by evaluating the distribution for different values of candidate "muon" masses. We simulate two 
GMSB-like models to show that this procedure can be used to detect the long-lived sleptons and 
measure both their mass and the neutralino mass. 
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I. INTRODUCTION 



Throughout the past four decades, collider experiments have unveiled the full structure 
of the Standard Model (SM), by gradually discovering a series of new, unstable particles 
through their decays to known light particles. It is perfectly plausible however that new 
particles beyond the SM are stable on collider-detector time-scales, by virtue of some new 
symmetry that forbids or suppresses their decays. If such new particles are charged, they 
must decay on cosmological time-scales, potentially leading to rich and testable dark matter 
scenarios Indeed, many extensions of the SM predict Long-Lived Charged Particles 
(LLCPs) that would not decay in the LHC detectors. These include for example sleptons in 
gauge-mediated or gravity-mediated supersymmetry-breaking models with a light gravitino, 
or KK-leptons in extra-dimension models (see for example the review [2J and references 
therein) . 

If produced at the LHC, an LLCP would traverse the entire detector, much like a muon, 
but unlike a muon, its speed would be smaller than one. There has been significant 
progress in recent years on (i) modifying trigger algorithms to ensure that LLCPs are not 
missed altogether if they are very slow, (ii) developing techniques for LLCP detection 
The latter, which have already led to the exclusion of various LLCP models jf!-^, mostly 
rely on the low LLCP speed in order to distinguish it from a muon. At the LHC, LLCPs 
could be produced either directly, or in cascade decays of heavier particles 1 . A large fraction 
of the latter may have close to one 2 . These fast sleptons will be hard to distinguish 
from muons based on Time-of-Flight methods. Fortunately, the fact that such fast LLCPs 
typically originate from the decay of a heavy particle, can provide new handles on them. 
Consider a long-lived slepton I which is produced in association with a lepton I from a 
neutralino decay, x° ~^ If the slepton is fast, it can be mistaken for a muon, so that its 
4-momentum is not measured correctly. The (mis)measured invariant mass of the opposite- 
sign (OS) lepton-slepton pair is therefore different from the neutralino mass. Still, if the fast 
slepton originated from a boosted neutralino, its mass is small compared to all the energy 
scales characterizing the decay, so that the mis-measured invariant mass distribution would 
peak somewhat below m x 2 — m 2 . This peak may be observable above the background. 
Furthermore, the peak can be enhanced by scanning over the possible slepton masses in 
order to obtain the correct neutralino peak. Taking all the OS "muon"-lepton pairs (where 
the "muon" could be either a real muon or a slepton LLCP), we can assign the "muon" 
a trial mass mscAN, and calculate the resulting "muon" -lepton invariant mass. Varying 
Tuscan, we will find that the highest peak in the "muon" -lepton distribution is obtained 
when mscAN coincides with the true slepton mass. If this peak is indeed observable over 
the background, one can isolate the LLCP candidates. In the process, both the slepton and 
neutralino masses are measured, the first from the best value of mscAN, and the second from 
the position of the invariant mass peak, calculated with the correct slepton mass. Finally, 
since at the end of the process the slepton and lepton 4-momenta are fully measured, one 
can hope to obtain information about the slepton and neutralino couplings and spin, as well 
as reconstruct the full event [9|. 



1 Thus for example, long-lived sleptons could come from Drell-Yan pair production, from cascade decays of 
strongly interacting particles, or (if the masses of these strongly interacting particles are very large), from 
the electroweak production of gauginos, followed by the decay of these gauginos to lighter sleptons. 

2 This would be even more of a problem if the LHC eventually goes to 14 TeV. 
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The applicability of our methods depends of course on the flavor of the LLCP. As a 
first step, we will study here a GMSB-like example with degenerate, long-lived right-handed 
sleptons, and consider only "muon" -electron pairs, which have little SM background in 
comparison with di-muon pairs. On the other hand, compared to models in which a single 
slepton is meta-stable, our model results in a factor of three reduction in the efficiency for 
sleptons. Furthermore, uncorrelated stau-electron and smuon-electron pairs contribute to 
the SUSY combinatorial background. It would be interesting to extend our approach to 
models with a single LLCP slepton, whether it is a stau, a selectron or a smuon, and we 
plan to investigate this in more detail in the future. 



II. THE SLEPTON-LEPTON MIS-MEASURED INVARIANT MASS DISTRIBU- 
TION 

When a high /3 slepton goes through the muon detectors, it can easily be mistaken for 
a muon. The slepton 3-momentum is then measured correctly, but the energy is computed 
from this 3-momentum assuming a zero slepton mass. Consider then such a mis-measured 
slepton, which originates from the neutralino decay x° ~~ ■ The measured, or rather, 
mis-measured, 11 invariant mass is, 

™% = (Pi + Pi) 2 , (!) 

where pi is the lepton 4-momentum and pj is the mis-measured slepton 4-momentum, ob- 
tained from its true 4-momentum by setting the slepton mass to zero. 

If the slepton is very fast, as would be the case if its mass mj is much smaller than the 
neutralino mass m x or the neutralino energy, one would expect the rh^ distribution to peak 
somewhat below the neutralino mass. Indeed, to leading order in 1 — /3 2 , where (3 is the 
slepton speed in units of c, takes the simple form, 



2 - " 2 -mf) (1- Tl *r^\ ) > ( 2 ) 
' m x + mf 1 — r cos 9 J 

mJ 2 — m? 

r = 2^ 2 3 

m x z + mf 

and 9 is the angle between the slepton and neutralino direction in the neutralino rest-frame. 
We see from fl2]), that m^ is a non-negative monotonically decreasing function of cos (9 whose 
maximum is obtained when the slepton and the neutralino are collinear, at m~. ~ m x 2 — mf 
in the fast slepton approximation. We can also invert this equation to obtain cos 9 in terms 
of rhjt, yielding the form of the m~ u distribution, 



m li = K m x 



- 2 1 + cos 9 



where 



dT 



mJ 2 + mj 2 ( 2m v 2 mr 2 
1 ± a — l - 11- 

m v — Tnr 




(4) 



where a < 1, and the sign is determined by the neutralino polarization (if the neutralino 
were a scalar, one would have a = 0). Thus, the rate is a monotonically increasing function 
of rhft, and peaks near its maximal value which is roughly m x — m 2 . In the Appendix, we 
show that this behavior persists beyond leading order in 1 — (3 2 , but the maximum allowed 
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value of m~ a is in general given by, 

r 2 ~\ 

(5) 



K max = K - <\ 



(E x + \p x \y 



which coincides with m x 2 — m- 2 in the limit of infinite neutralino energy 3 . 

We can thus try to use the m~ n distribution to detect the slepton LLCPs. The OS 
"muon" -lepton invariant mass distribution behaves differently for "muons" which are true 
muons, and for "muons" which are meta-stable sleptons. The latter have an m- u distribution 
which is monotonically increasing and exhibits a peak at the maximal allowed value of m^. 
In contrast, the invariant mass distribution of real muon-lepton pairs usually falls off for 
large values of the invariant mass, unless of course the muon-lepton pair is correlated, and 
comes from the decay of some heavy particle, such as the Z. Such backgrounds are however 
reducible, and can be eliminated, unless of course m x 2 — mf happens to lie close to the Z 
mass. 



III. THE m n DISTRIBUTION IN A GMSB-LIKE MODEL AT A 7 TEV LHC 

In order to demonstrate our techniques, we simulate two GMSB-like models for a 7 TeV 
LHC with the SUSY-breaking parameters 

A = 5 x 10 4 GeV , M msg = 2.5 x 10 5 GeV , N 5 = 3 , tan/3 = 5 , sign/x = + , C grav = 5000 . 

(6) 



Here A is the ratio of the GMSB-singlet F term to the messenger scale [lOj, llll], N 5 is 
the number of messenger multiplets, and the parameter C grav is the ratio of the gravitino 
mass to what it would be if the only source of supersymmetry breaking were the GMSB 
singlet F-term [12|. Since the latter F-term is typically generated from larger F-terms, 
such large values of C gmv are quite plausible, resulting in heavier gravitinos compared to the 
naive GMSB estimate, and therefore in long lifetimes of the Next-to-Lightest-SuperPartner 
(NLSP). For this choice of GMSB parameters, the NLSP is a right-handed slepton 4 . We list 



the superpartner masses in Table fl] These masses are calculated using SPICE [13J , which 



is based on SoftSUSY [14| and SUSYHIT [15j. Since we are interested in fast LLCPs, we 



modify this model by taking the three right-handed sleptons to be degenerate with a mass 
of either 110 GeV or 130 GeV. We will refer to the two models as the 130 and 110 models, 
according to the NLSP mass. 

We generate a total of 10000 SUSY events, corresponding to an integrated luminosity 
of 146 fb _1 at a 7 TeV LHC. These come from strong-strong, strong-weak, and weak-weak 
pair production, which at leading order have a cross-section of 68 fb. We do not include 
Drell-Yan slepton pair production, which has a comparable cross-section, since most of the 
resulting sleptons do not pass our selection cuts. 

The dominant relevant SM background comes from ti for which we have generated a 



3 Note that rhr, > follows from the fact that pi and pj are both lightlike. m ; ~ ; = 0, however, can be 
obtained only for sufficiently energetic neutralinos, with E x — \p x \ < m~ v 

4 Recall that the gaugino masses grow as N$ , while the scalar masses only as \fN§. Already for N$ = 3 the 
neutralino is heavier than the right handed sleptons. 
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TABLE I: Spectrum of the GMSB model calculated from SPICE. The three lightest slepton masses 
are set by hand to 130 or 110 GeV. 

sample of 420000 events corresponding to a luminosity of 5 fb _1 at a 7 TeV LHC. The ti 
sample is produced at leading order with a cross-section of 85 pb. Unless stated other- 
wise, all the results we show below are normalized to 5 fb _1 . To generate events we used 
MadGraph-MadEvent (MGME) [l^|, with FeynRules [p7j ]. The resulting events are decayed 
using BRIDGE and put back into MGME's Pythia-PGS package [342l| which includes 



hadronization and initial and final state radiation. The slepton and lepton momenta are 
then smeared according to reported experimental resolution 5 . 

In Fig. [T] we show the /3-distributions of the sleptons in both models, taken from truth 
(i.e., before any smearing). Indeed, most of the sleptons are fast. For both models, roughly 
50% of the sleptons have j5 > 0.9 and 25% have (3 > 0.95. As explained in the Introduction, 
we are interested in fast selectrons which are mistaken for muons. We propose to detect 
these selectrons through the selectron- lepton mis- measured invariant-mass rn~ u . The m~ u 
distribution of signal events (i.e., OS ee pairs coming from the decay x° ~^ e ± e =F ) is shown 
in Fig. [2] for the two models, with a pt cut of 20 GeV on both the electron and the selectron. 
As expected, the distributions increase monotonically and peak near the maximal allowed 
value of m| ~ m 2 - m- 2 (see ©), which is about 150 GeV (165 GeV) for the 130 (110) 
model. 

The relevant backgrounds include uncorrelated OS slepton-electron pairs from SUSY 
events, with the sleptons in this case being either selectrons, smuons, or staus, and OS 



5 Note that we do not use the PGS detector simulation. 
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(a) Normalized /3-spectrum of sleptons for (b) Normalized /3-spectrum of sleptons for 

m r = 110 GeV mj = 130 GeV 



FIG. 1: The true beta spectrum (unsmeared momenta) of sleptons. 




(a) m f = 110 GeV (b) m ; ~ = 130 GeV 

FIG. 2: The frige distribution of OS ee pairs from the decay x° ~~ > e^e^, with pr > 20 GeV for 
both electron and selectron, for 146 fb _1 , for the two models. 

muon-electron pairs from both SUSY events and SM production. The main source of the 
latter is tt production. To enhance the signal over the background we impose the following 
cuts: 

1. pt > 20 GeV on both "muons" and electrons. 

2. p T > 100 GeV and p > 250 GeV on the "muon" . 

3. cosOjt > 0.6 where 9~ u is the angle between the "muon" and the electron. 

Here and in the following, "muon" denotes either a real muon or a slepton LLCP. The number 
of events surviving these cuts in each of the models and in the tt sample at 5 fb _1 are collected 
in Table [Til The combined set of cuts significantly increases the signal-to-background ratio. 
Still, for the models we study with 5 fb _1 , the signal peak is not big enough to dominate the 
background. In Fig. [3] we show the OS e/x invariant mass distributions: signal, background 
(SUSY background and tt), and signal plus background after the three cuts for the two 
models, for an integrated luminosity of 5 fb _1 . In the signal distribution "muons" are 
selectrons, whereas in the two other distributions "muons" are either real muons, or any 
metastable-slepton (selectron, smuon or stau). Although the signal+background distribution 
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TABLE II: The number of events surviving the set of cuts described in the text for mr = 
110, 130 GeV and ti normalized to 5 fb _1 
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FIG. 3: The invariant mass distributions of OS e\x pairs after the three cuts in the m? =110 GeV 



(left) and mr = 130 GeV (right) models at 5 fb~ . "Muons" in the signal distribution (red, dash) 
are selectrons. "Muons" in the background (black, dash-dot) and signal+background (blue, solid) 
distributions are either muons or any metastable-slepton (selectrons, smuon or stau). 



exhibits a small peak near m x 2 — mf, it does not clearly indicate the presence of a decaying 
particle. In the next section, we will see how this peak can be enhanced. 



IV. SCANNING FOR THE SLEPTON MASS 



As we saw above, even after our event selection, the signal peak is small compared to 
the efi background coming mainly from ti production. It is clear however how to enhance 
the peak: if we evaluate each "muon" four-momentum with the correct slepton mass, the 
signal distribution would peak at the neutralino mass, while uncorrelated fi — e pairs 
would not. As a result, the signal peak may be visible over the background. We therefore 
study the electron- "muon" invariant mass distributions adding a trial mass tj-scan for each 
candidate "muon". For each mscAN, we calculate the maximal bin height, CmscAN- We 
expect that the maximal Cm SCAN wou ld occur when mscAN coincides with the true slepton 
mass, since then the signal should exhibit the narrowest peak, centered at the neutralino 
mass. In Fig. H]we show Cm SCA N as a function of mscAN, varying the latter over 300 values 
in the range 100-400 GeV, for the dataset obtained for the 110 and 130 models. Denoting 
the value of Tuscan at which Cm SCA N i s maximal by ^s C an; we ^ n< ^ ^ na ^ ^ n eacn °^ ^ ne ^ wo 
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(a) mj = 110 GeV 



(b) rrij = 130 GcV 



FIG. 4: The £ measures for each sample mass /n-gcAN) for the two models. 



models, Tuscan precisely reproduces the true slepton mass: rnl CAN = m e- We can now use 
m scAN t° plot the "muon" -electron invariant mass distribution, calculated with mg CAN for 
each "muon" candidate. The results are shown in Fig. The signal peak is clearly visible 
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FIG. 5: OS "//'— e Invariant mass distribution constructed with the for each candidate 

"muon" . Events include both signal and background after cuts. 



now over the background. 

There are a few things we can learn from these plots. First, the existence of a massive 
LLCP can be established, and its mass measured to be ^scan- Second, one can now isolate 
a signal-rich data sample, which includes "muon" -electron pairs whose mg CAN -corrected 
invariant mass lies close to the peak. Each "muon" in this sample can be re-examined, 
to see whether its track can be distinguished from a true muon. Third, we learn that the 
LLCPs in the signal originate from the decay of a heavier particle. The mass of this particle 
is given by the location of the peak in Fig. Finally, with this information, the signal 
events can be fully reconstructed, allowing for measurements of the masses and couplings of 
heavier particles with selectron and neutralino daughters. 

One could worry that by evaluating the distributions with some arbitrary value of it-scan j 
combined with the cuts we use, we are artificially introducing some features into the dis- 
tributions. Thus for example, for uncorrelated OS muon pairs, with a cut of p^ nt {pl nt ) on 
the momentum of the slower (faster) muon, one would expect the majority of the leptons 
to have momenta close to the cut so that the peak of the mgcAN-corrected distribution is 
obtained for 

m?r scAN - m\ 2 ~ ml cm (l + p^/p^) , (7) 
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where m\ 2 is the true invariant mass, and m('™ SCAN is the mgcAN-corrected invariant mass, 
evaluated with a mass m SC AN for the faster muon. 

To address this issue, we take two different approaches. First, we examine the behavior 
of the distributions for both wrong values of mscAN and for mscAN = m scAN- This is done 
in Figs. [6] and [7] where we plot the distributions for 16 values of mscAN for each model. The 
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FIG. 6: Invariant mass distributions for each sample mass with mj = 110 GeV 

results suggest a two-bump structure, with the lower bump presumably corresponding to the 
background and lying somewhat above mscAN, consistent with ((7|). However, for values of 
"^scan that are very different from the true slepton mass no distinct features are observed. 
Second, we examine the "raw" and mg CAN -corrected distributions of the SUSY back- 
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FIG. 7: Invariant mass distributions for each sample mass with mr = 130 GeV 

ground and of the ti background separately from the signal (see Fig. |S]) 6 . Indeed, it is easy 
to see that including mg CAN does not generate any new features in the distributions. In 
fact, the one feature that does appear in the SUSY background upon including mg CAN , is a 
peak at around 350 GeV, which is nothing but the Due to the low SUSY cross-section, 
however, this "peak" contains only one event.. 

Note that the cross-checks of Fig. [61 Fig. [?l Fig. [8] can be done in a real experiment. 
This is clearly true for the first two, but also for Fig. [HJ To that end, one can isolate the 



6 We only show the results for the 110 model. The results for the 130 model are very similar. 
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(a) SUSY background with massless "muon" (b) SUSY background with massive "muon" 

candidates candidates of mass rascAN = 110 GeV 




,M e(I ,M e|I 

(c) ti with massless muons (d) ti with massive muons of mass 

toscan = 110 GeV 



FIG. 8: The OS e/i invariant mass distribution of SUSY background (top) in the 110 model, and 
of the ti background (bottom), after the three cuts. The left panels show the raw distributions, 
and the right panels show the distributions evaluated with mscAN = 110 GeV. "muons" include 
both sleptons and real muons in the SUSY background and only real muons in the ti background. 



signal-rich sample and the background-rich sample as described above, and plot the raw and 
mg CAN -corrected distribution for each one of them separately. 



V. CONCLUSIONS 



We have studied a method for detecting fast long-lived sleptons and for measuring their 
mass. This method relies on the mis-measured invariant-mass distribution of slepton-lepton 
pairs with the slepton faking a muon. As a first step, we have only studied here a model with 
three degenerate metastable sleptons. This has both advantages and disadvantages. On the 
positive side, the model has a metastable selectron, which results in selectron-electron pairs 
produced in neutralino decays. These produce observable features in the "muon" -electron 
channel which has very low backgrounds from correlated lepton pairs. On the down side, by 
using the "muon" -electron channel only, we are keeping only a third of the signal, compared 
to models with a single metastable slepton. Furthermore, the metastable stau and smuon 
can also be mistaken for muons, and therefore contribute to the background. 
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The models in this work were chosen especially for their low cross sections, so that they 
are not already excluded by time-of-flight or specific ionization methods Though our 

analysis is only a "leading-order" analysis, it is tempting to estimate the lower limit on the 
cross-sections of detectable models. From Fig. it is clear that much lower neutralino mass 
peaks (by approximately a factor of 5) could still be discriminated from the background. 
The efficiency of the method, however, depends on the neutralino momentum distribution 
which in turn is affected by the mass spectrum, and is thus model dependent. 

The methods we describe would become even more important in a 14 TeV LHC, since 
then the slepton /3-spectrum will peak more sharply near /3 ~ 1. 
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Appendix A: Mis- measured Invariant Mass Distribution in a 2-Body Decay 

In this Appendix, we will derive the functional form of the m~ n distribution, and show that 
it is a monotonically increasing function, and therefore peaks at the maximal value of m^. 
Note that, although the mis- measured slepton momentum p^ is not a physical momentum, 
it is nonetheless a 4- vector so that m~ (see (JT])) is Lorentz invariant and the neutralino 
differential decay rate dT/drh^ can be calculated in any frame. Although it is natural to 
think of rhfi in the Lab-frame, since this is where the slepton 3-momentum is measured, it 
is instructive to express it in terms of the neutralino energy E x (in the lab-frame) and the 
angle 9 between the slepton and neutralino direction in the neutralino rest frame, since this 
angle does not depend on the neutralino energy. 

We consider the decay x° ~>* ■ m the neutralino rest frame, the particle momenta 
are: 

Purest = {Ej Rest , -Pi Rest) Pflrest = {pi R est ,Pl Rest) (Al) 

where p denotes the slepton momentum, and, using energy-momentum conservation, 

m 2 x — mf m 2 x + m 2 

Pl Rest = 2m x El Rest = 2m x ' {M) 

Taking the neutralino direction to be +z, and writing the lepton and the mis-measured 
slepton 4-momenta in the Lab-frame as, 

Pl Lab = (PhPl Lab) , P 11 Lab = (P,P Lab) (A3) 

where 

P = |PLab| Pl = \PI Lab| , (A4) 
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we have 



Pi 



Pi Rest 



( E x+Px cos9 ) ■ 



There are several useful ways to write the mis-measured invariant mass-squared, 



™1 = (PiLab + P" Lab)' = (Pi + PY 



Py 



M-m 



2pi(E-p) 



(A5) 



(A6) 



where E is the slepton energy, E = \/p 2 + mj 2 . Note that m~ = (pi + p) 2 — p 2 x is a function 
of p x and cos 9. 

In the neutralino rest frame, we have 



dT 



2 2 

m x — m 



dcos8 32irml 

A. 



(A7) 



with |M| 2 = A(l ± acos^), where a is a positive constant with a < 1, and the sign depends 
on the neutralino polarization. 

We can write the differential decay rate in terms of rh~ lv using Eqs. (IA5[) and (IA6[) . 



dT 

difiTi 



mr 



^m x p x lrhl + \p x \ 2 



(E x - Jm\ + \p x 



To determine the sign of d 2 Y '/dm^ it is convenient to start from expression (1A7[) . and use 
the chain rule to convert derivatives with respect to cos 8 into derivatives with respect to 
m^. Doing that, we find that dT/dm^ is a monotonic function of m^. 



J. L. Feng and J. Kumar, "The WIMPless Miracle: Dark-Matter Particles without 
Weak-Scale Masses or Weak Interactions," \Phys.Rev.Lett. 101 (2008) 231301] 
|arXiv: 0803. 4196 [hep-ph][ 

M. Fairbairn, A. Kraan, D. Milstead, T. Sjostrand, P. Z. Skands, et al, "Stable massive 



particles at colliders," \Phys.Rept. 43 8 (2007) l-63[ |arXiv:hep-ph/0611040 [hep-ph] 



A. Nisati, S. Petrarca, and G. Salvini, "On the possible detection of massive stable exotic 
particles at the LHC," \Mod. Phys. Lett. A12 (1997) 2213-2222} |arXiv:hep-ph/9707376[ 
CDF Collaboration, A. Connolly, "Search for longlived charged massive particles at CDF," 
|arXiv:hep-ex/9904010 



S. Tarem, S. Bressler, H. Nomoto, and A. Di Mattia, "Trigger and reconstruction for heavy 
long-lived charged particles with the ATLAS detector," \Eur. Phys. J. C62 (2009) 281-292 [ 
ATLAS Collaboration Collaboration, G. Aad et al, "Search for stable hadronising 
squarks and gluinos with the ATLAS experiment at the LHC," 



Phys. Lett. B701 (2011) 1-19, arXiv: 1103. 1984 [hep-ex] 



ATLAS Collaboration Collaboration, G. Aad et al, "Search for Heavy Long-Lived 
Charged Particles with the ATLAS detector in pp collisions at sqrt(s) = 7 TeV," 
\Phys.Lett. B703 (2011) 428-446[ |arXiv : 1106.4495 [hep-ex]| 

CMS Collaboration Collaboration, V. Khachatryan et al, "Search for Heavy Stable 



Charged Particles in pp collisions at sqrt(s)=7 TeV," JHEP 1103 (2011) 024 



13 



|arXiv: 1101. 1645 [hep-ex] 



J. L. Feng, S. T. French, I. Galon, C. G. Lester, Y. Nir, et al, "Measuring Slepton Masses 

and Mixings at the LHC," | JHEP 1001 (2010) 0471 |arXiv:0910.1618 [hep-ph] [ 

M. Dine, A. E. Nelson, Y. Nir, and Y. Shirman, "New tools for low-energy dynamical 



super symmetry breaking," Phys.Rev. D53 (1996) 2658-2669 



arXiv:hep-ph/9507378 [hep-ph] 



M. Dine, A. E. Nelson, and Y. Shirman, "Low-energy dynamical supersymmetry breaking 
simplified," \Phys.Rev. D51 (1995) 1362-1370[ |arXiv : hep-ph/9408384 [hep-ph]] 
"ATLAS: Detector and physics performance technical design report. Volume 2,". 
G. Engelhard, J. L. Feng, I. Galon, D. Sanford, and F. Yu, "SPICE: Simulation Package for 



Including Flavor in Collider Events," Comput.Ph ys.Commun. 18 1 (2010) 213 -226 
|arXiv: 0904. 1415 [hep-ph] [ 

B. Allanach, "SOFTSUSY: a program for calculating supersymmetric spectra," 
\Comput.Phys.Commun. 143 (2002) 305-331] |arXiv:hep-ph/0104145 [hep-ph]] 
A. Djouadi, M. Muhlleitner, and M. Spira, "Decays of supersymmetric particles: The 
Program SUSY-HIT (SUspect-SdecaY-Hdecay-InTerface)," Acta Phys.Polon. B38 (2007) 
635-644, |arXiv:hep-ph/0609292 [hep-ph]] 

J. Alwall, P. Demin, S. de Visscher, R. Frederix, M. Herquet, et al., "MadGraph/MadEvent 
v4: The New Web Generation," | | JHEP 0709 (2007) 028] |arXiv: 0706. 2334 [hep-ph] | 
N. D. Christensen and C. Duhr, "FeynRules - Feynman rules made easy," 
\Comput.Phys.Commun. 180 (2009) 1614-1641] |arXiv: 0806 .4194 [hep-ph] 



P. Meade and M. Reece, "BRIDGE: Branching ratio inquiry / decay generated events," 
|arXiv:hep-ph/0703031 [hep-ph] | 

J. A. et. al, "MadGraph - Pythia Interface," http://madgraph.hep.uiuc.edu/ . 
T. Sjostrand, S. Mrenna, and P. Z. Skands, "PYTHIA 6.4 Physics and Manual," 
\JHEP 0605 (2006) 026] |arXiv:hep-ph/0603175 [hep-ph] 



[21] J. Conway, "PGS - Preaty Good Simulation ," 

http://www. physics. ucdavis. edu/ conway/research/software/pgs/pgs4-general.htm . 



14 



